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ie i &itmmm±.<Dm%tF ^^m-r^mm^mz^x, 

*<kx<Dmt LxmuMB-f z>m<7)tp<k<7>mf zw-h, 
kx^>x, mrnx^t h i>kKnmm<DUftmM±x<Dmmm<7)m*7E.&2>z- 

umm 3 j 2 &7tmmM±.<mm $ 2 ^tc«^^ ^ , 2 

IE 2 ^TC^fl it± <7) M £fc F £ TEirt £ fit HQ # V* K *3 v » X , 

<D2fe<7)m* J tti J ?tiD f D f ffi, DfS, Df^tU buIED f m<bHUfE 

ffi^u HuiED f ffitfrisD f it&mfftom&iwu^ mn-%v>Wi&i**<o& 

*Mjt£>/h£v>*>0£D tmitt U li>rfED f mffitmifED f m^bOM^$r*LHfr 
IEJ&¥ <7> K £ v> xm*¥ <7> $*<kX<D HtffE^^ f <r>m £ #o^p® LX^rti 



mwEW2 003-3067030 
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3uao^-co f <om*Tmb%^£ ■? tcffisD f hcm^'buIbd 1 jy. 
<Dft®m%-m<v2fe<Dm* j tnmi) 1 1, Dfi, Df*, Df^L Df^, 

(4-trn, [WJ^^^{4^0^tH<7)/h^v^<7)^D f ffizjfci: U fufeDf^h 
f^mt U tuIED f mStHtlfBD f itthHyfBD f § Lliufejfe¥ 

tuf EJ& ¥ o^km^[6] O^T- * K £ it * ¥® £<7>{trt ? lfe ¥ o *<fr<£Hg i^i 
^7g^-COff^x.^v^^;, HufBDfmM, HfrfBD f IfrlfcWfT 
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ED f 1 &TO**:tf>56»fcJ: OfcBrtfrtS-fcfcJ: ^IL^IU 

^bH^ifijo^T-^, mnttfa<Dfe*g„ mytm-xmote*^ At 
m&^ft^t&i 1 *. imtoiji*wtifr&„ mnmii^m^^ sow^i^i 

w-m^mm^x^mmmmmmt Lt^^fMttai^Mo 

tf^-COffifc LTHtrfE^Bf &fl&*>*'k«H*f *«f*>. il^^oiw 

mum 7 3 2 &7tnm&.-te>m& $ntz2 xjtm&fcf- 1 . m 2 $otm 

ffirlBMSfef OM-«3l^ll«)2fSO|tS:^-e*nDf 3K, DfE Dfi, D f 4b 
£:U MED f miifffBD f ISibWW-Znm&ltW, mW^<r>%r&l±*<7>1&tt 
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3m<D^^X(DMUm%kf <DW.ZT®b%:^*£ 7 HlfffED f ICMWliiffED ft 

Ft Lxfc&z^&tzmz-tz*z&m-rznftmwio 

itmm 8 3 3 <k?tmmm±v>m& s ntz 3 $otmm.&i- 1 , m 3 ^« 
^te^^us ^ii* is v > x m t & m k * o x £ n * m 3 ^7c«it^T-± 

IE 3 &7tmmm±<o m $k f * ^m-r z mm * mft-r %> mnmm k & ^ r , 

<Dfr®mfrm<0 2fe<7>m*-Ztl J ?fftD f I, DfB, Dfft, Dft, Df^, 
DfitU frfBD f y&tmUD f S^flF-^Ofc-frli-tfo, HffiHf-4>*£-«:* 
<D&j$m<7)s\s2\,*i> <D&D imst U mflSD f mtmiED f 
li-tfn, ra^O*^(i-€-0^1it<^/h$v^o*D f mJbfc L, HtrfBD f ^ <h 
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^7l&Wt'it^ f d BuffiD f luffiD f ^bS.^Btf 

BSD f lJ^TOf^^tia »9^^-t-^C^tcJ; ^TO^£^H-f£ 

n<v7m<7>*<L^<Di mzTm^z^x^^mmD f mn. frieD f mtwm 

imttfovfe^g^ mmttfa<Dfc*&. si^isi^^, ft 
m*^fS]<7)^T-*. itwmjjfa<7>fc=?-&. mnm^<D^^ Rimyn-m^ 

is i ^7cnnit± o m$t f * & m r% * & ^ n > m ^ - * -e 
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umm 103 2 >xjtnmm±<Dffij& z ntz 2 >kitmm^ m 2 
mm. fe? <?)&m *<c^ & v > x z & * k x o r $ *i & ^ 2 ^7c#is»r- 

lWTSB 2 &x*3Wft±tf> 08 F £ fcM-r Z> fit Wite * 1" * # n > kf a - * 
<7>2fe<7)mi*ti*rtLD f DflS, Dfg, DfttfcU luffiD f mi:fff£ 

d f istimw^m&infu, m^^m&it^^mttm^^z^^nzv f ^ 

*MItf>'h$^fc<0£D f flMbfc U huIBD f t mrfSD f mt<DM£ Lm 

§^»£-«±, Hul£¥®^<7)^bm^[S]^^T*-e<7)ji^\ iM&^£&#*rl-*ife 
3jfto*Nye<*> f (DrntTMhtc^x 0 izmuD fmmRxrmmD mat* 1 \>x 

[if « 113 3 ^HS«±OiS $ *Lfc 3 ^7C«^f&T- , MS>kjt 
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0>tllll^i0 2f&0|fiS:«etL-r*LD f JPu D f H, DfS, Dft> Df^ 

Dfitu ttriBD f m^ButeD f s^m^^^ti-tfn, m^<vm<£tet 

(i-tfn, liJIf-i-o^^a-eo^ll^/h^v^o^D f mitt U mmDi^t 
iJO&t L, tirfSD f mM^HUfSD f m^^HufED f ^oM§**Lmrie^¥ 

IBD f 1 &T<7)m±<7)7£.mz£ ^^.Wc^t^^tizX ^TO^flL, 

^<7)7^4 1 ^-e<7) f lit^Tll'b^v^J;^ tcplED fIS, lufSD f M^mr 

ifiD f i&T<7)m±<D7£mz± ^mirf tic* ^¥®z:£^hu 

fc*ffc*. m^^o^*, iro^ftof&T-*, mm^ft^*&?-,£u it 
m%^ftof&^au itM^ft^?-*, itM^ft^^^ &£mm*fc;w 

[o o o 1 3 
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[0 0 0 2] 

[0 0 0 3] 
[0 0 0 43 

[0 0 0 53 
[0 0 0 63 

£'j -^><D&%rp<D]&.W)Jjfe£:<omJ: TUeber die Fortpf lanzung ebener Luf 
twellen von endlicher Schwingungsweite (Abhandlungen der Koniglichen Ge 
sellschaft der Wissenschaf ten zu Gottingen, 8 (1860), 43-65) J tcH/^-£\ 

[0 0 0 73 
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iO'J-7^W«©fti»ffitLT, (Total Valuation Di 

minishing&, ^TTVDffitlB-t) tPf tfn*^ffi«±^b^^- A, fUMBaifc 
SrSAi-*^KW^NB»ffi (Essentially Non-Oscillatory J^TEN0y££ l£-f 
) ^^-^-P-h'/jfL'^ (Monetone Upstream Scheme for Conservation Laws^ lck~f 
MUSCL&fcK1-) . lEfitfMffe (Positively Preserving, &T??&tttt)m* <D 
Rt^^ft^SL, -€-*Lt>ti»aiT*ofe«9, il6JRt*o/:i)LO^o 
[0 0 0 8] 

a F(x,t) / at =ac(x, t) f(x, t) / ax 

... (D 

#x.&o F(X, T) t LT, X.Tt±^T«ft«F^n^jS»tmfkl*W*»-t<>Ofc 

-r 

[0 0 0 93 

5£ (1) *i<IBI«WL^ ±^Lf:i^:, tfit*ffcfc¥K:i3V*T«:ft* 
fctf-JMH* «S $ tLT 0 TVD& t it > TVD&fr 

TV(F[ • ,n])) = 2j IF[j+l,n] - F[j,n]| 

- (2) 

TV(F[ • ,n+U) ^ TV(F[ • ,n]) 

- (3) 

[00103 

j£^^G£T<7)'t£lf#P±»fLv£ (Monetone Upstream Scheme for Conservation 

Laws , MuscuitfBi-) (i. r- ? omK&t&v) —?>mmt<D 2 ^frbffi- 

[00113 

MUSCLE t TVD& £ £ d <h *fc M»tc^a«fc*« CWSPfc-faSteU 
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[0 0 12] 

[0 0 133 
[0 0 143 

H[I,J](x) = F(IJ) + G(I,J)x 

... (4) 

^ftitv^t Lrm^mm^M^m^m^^o i> l*>m&c(i+i/2,j) 

Q(I+1/2,J) := C(I+l/2,J)H[IJ](Ax/2 - C(IJ) At/2) 

- (5) 

t u ^-e&u-tui\ 

Q(I+1/2J) := C(I+l/2,J)H[I + l,J](-Ax/2 - C(I+1J) At/2) 

- (6) 

F(I,J+1) = F(I,J) + At (QCI+1/2, J) - QCI-1/2, J) )/A x 

- (7) 

ti-zztizx f), ^ (i) *miwi<zM< :u j t^o 

[0 0 153 

flidj. ^ (4) OE^l-^Mm^i3V^T^lSe^^-r^n-yMi$:G(IJ)(im 
[0 0 163 

G(I,J) = AF(I,J)/Ax 

- (8) 
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[0 0 173 

Sp( I, J): =F(I+1,J) — F(I,J) 
Sm( I, J): =F(I,J) — F(I-l.J) 

- (9) 

*m\-f%>o l^CXTVDWIRHSfct LT. Minmodt D^tf *L& & ©|:i5V»tli, AF 
(I,J)t±, 

AF(I,J) = (sgn(Sp(I,J) + sgn(Sm(I,J)) )Min(|Sp(IJ) I, ISm(I, J) I ) / 2 

(10) 

AF(I,J) = (sgn(Sp(I,J) + sgn(Sm(I,J)) ) 

Min((ISp(I,J)+Sm(I,J)|/2 ,2|Sp(I, J) I, 2ISm(IJ)l ) /2 

-(11) 

•C * 0 > Superbee t P^lffi & c7) t? (± , 

AF (I,J) = (sgn(Sp(I,J) + sgn(Sm(I,J)) ) 

Min(Max(|Sp(I,J)|, ISm(IJ)l) ,2|Sp(I,J)l, 2|Sm(I,J)| ) /2 

- (12) 

t*4Hm„ lit. sgn o liftias-e^-So -eo^i^iE 

[0 0 18] 

C tl <b (O Z t (i , TA. Suresh, "Positivity-Perserving Schemes in 

Multidimensions" SIAM J. Sci. Comput. 22, 1184-1198 (2000)J iZUHBL^fiX 

[0 0 193 

^IH7~I1I9 i:*^tv>5 c *7^7iU 0 0 00^0 

Sfltfc, Wm^77C|L/;o HI 7 35«Minmod, IH 8 **¥*&IWIRM», H9 
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# ? Superbee £ Vfr l£it & & <D X $> Z> 0 ^ftfbOEI^^^&J:?^ t&Mfctkff 
[0 0 2 0] 

[0 0 2 1] 

[0 0 2 21 

Vmin = Min( F(I-1,J+1. T)-F(I, J, T), F(I,J+1, T)-F(I, J, T), 
FU+l.J+l, T)-F(I,J, T),F(I-1J, T)-F(I, J, T), 

F(I+1.J, T)-F(I, J, T), 
F(I-1,J-1, T)-F(I, J, T), F(I,J-1, T)-F(I, J, T), 
Fd+lJ-l, T)-F(I, J, T) ) 

Vmax = Max( F(I-1,J+1, T)-F(I, J, T), F(I,J+1, T)-F(I, J, T), 
F(I+1,J+1, T)-F(I, J, T), F(I-1,J, T)-F(I, J, T), 
e, F(I+1,J, T)-F(I, J, T), 

F(I-1,J-1, T)-F(I, J, T), FCIJ-I, T)-F(I, J, T), 
F(I+1,J-1, T)-F(I, J, T) ) 

- (13) 

Sx(IJ, T) = (F(I+1J, T)-F(I-1,J, T) )/2 
Sy(I,J, T) = (F(I,J+1, T)-F(I,J-1, T) )/2 

••• (14) 

[0 0 2 33 

V = 2 ( MindVminl, I Vmax I) 

/(ISxd.J, T)| +ISx(I,J, T)l+e) 
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••• (15) 

Gx(I,J, T) = Min(l,V)Sx(I,J, T)/Ax 
Gy(I,J, T) = Min(l,V)Sy(I,J, T)/ Ay 

••• (16) 

tL^i:, (I, J) »^F«30«mM»H(I,J,T)(x,y)S: 

H(I,J,T)(x,y) = F(I,J, T) + Gx(IJ, T) x + Gy(I,J, T)y 

- (17) 

[0 0 2 4] 

dftfcfijffl-r&fc, no (a) izTF-tmmm^ttLx. mio (b) n^-r 

[0 0 2 5] 

mwtm®: L£d t-tz, mm] 

[0 0 2 6] 

2ffiO*JE»ttSlt#Ui5tt4at##«i* (VOFft) &c£v>T, ff^SrMTK 
$ tz t i m ffl^m H4ji>kv»9|l(IItt, § * IH * k *> 

ra^-r. ^*o^-e(±. mzfm&ffiiix%tzh frizz 

oT$> < Mfilfrfr-ofzo 
[0 0 2 73 

mx*)*), ^(DtztbK, — ^TCH^-r^T^^UXAtcov^TJi, ±ifi:£ (4) 

^Til'^'JXACov^li, ±IB^: (17) K;J3lt&1$WM£c&<fc 9&^*>'<*>&;: 
M # £ & - £ ^ £> o £ 0 
[0 0 2 83 

ttilE# 2003-3067030 



#SS 2002-223145 



-^-vl 14/ 



[0 0 2 93 

*^<toT^^tt^>iti^7c«^T-±oMmf ^jams*^^ mwm 

[0 0 3 0] 

^^(7)0- <7)ffip B ^ritrf±> 2 ^7E*fItiSUi<oM;£ fitz 2 ^7n«5t^t 
iiKJ:^ BtjfB2^7cllfIlt±OP^i5:F^^«-t^ffiP^^tlioV^T, mB-?2> 

^\ ^T-***^r1--i>fla¥^i-<7)3J3a^4 , >L^-co f <7)ft£j@;i&^J; 9 IZD f 
mB&tfD f aMt* 1 WTOf^O^tCi «9 5e»fg-r*-fc^«t 0 iftfe^ffiZ. 
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It* 1 TOi^tCJ; *) fc&fe-f & Z. t iz £ Ofuffi^MZl^mU #<b 
[0 0 3 13 

> it 3 ^Ttmm^^^m^^z^xm^i zmizx-o x^m^^mz >x 

<Dmum$kf <vK®mftm.<r)2fe(Dm.* j tti : zffr'D f m> Dim, Dfs, Df 
it, Df^ Dfftti, tuiBD f mtMiaD f m-timw^com&izifv, mw 
■%-cr>m&it j tvmm<7)'bz^i><DZD fistu mmv f m^tuieD f 

tfffSD f ^iifjfBD f iWlifr^^ii^n, ra^oi#^(i-e<0^t'fit(7)/h 
£^&<7)£D f U mifBD f mM^HflfSD f ffitt t buIBD f »II 

itRi/muD f i\>XT<Dm-K<Di£.mz& VTzmm-rzztbzx v^wz,* 
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rtToTt^^, m^t^o^*, mmttfte>&F&, mm^fsjo^ 
^m^ft^^*, ^HHfe^ffi]ot§T-^> nmm^<D^,^ wn^. 

[0 0 3 2] 

£K«fc 0> mfBl^x^IM±ora^:F^^i-^iir B Wvi*^T-r^ll^S 

& & £ £ ic im § £ -if n K&fe l , i£*&T-_t o * * mmmft t mmm^w^ 

[0 0 3 3] 

*mw<7>mr.<7)mnmw:'?iz, 2>X7tnmM±<7)mj&zntz 2>kjtmm^t 

^ it2^x«^^0#|&4 , 'L4;i3V^T»^f ^^tcj:or^a^tL^M2^ 

, Df^tU HuieD f m^HuffiD f ffi^ll^O^fi-trn, [Btff-^O 
>eo*fe*HI<7)/h£v^(7)£D f Mnt U HtjfSD f mtHyfSD f Ib^H^O^ 
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[0 0 3 43 

> m 3 ^xf#^£^#flS^L^£^Tlg^2 X o r^a £ n£fg 3 >k 

rti&fr^tzmiz^ MwmtJiMJj&iizm&o 6!&^^lt. huiemo^l^ 

t^6^^>L^i:<7)BufSMmf OM-fflM^fll^2^0flI^-eti-etiD fl, D 
fR Dfl%, DflL Df^ DflH, buIBD f m^BUlSD f 

f ffii:HtrfSD f »^^0^(i-tfn, |5]^co^(i-eo^ffl^/h$v^ 
<DltD f Wilkb L> mflSD f ^£ mifSD f ife^llft^O^te-lfn, mW^^m 

^ii^<Dmm^^^^i>^iD i^.mt u huied f me^tusED f mttm 

BED f »Il^f LMMm¥<7)*^!<z&l}Z>m%k f ^tit^ii&M&gifc^fgl 

ttilEtf 2 003-3067030 
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d f mm. m?>D f mtRz/muD fjzmzi ^Tot^^^si: «t •) ^^i- 

nmmD f misD f m^o^isd f 1 i^r<DM±^mzx ose 
m^mofe?*, wm^ft^iti^ imm^m^m^^ im^n^^ 

[0 0 3 5] 
[0 0 3 63 



tBIE#2 003-3067030 



#0 2002-223145 



^-v: 19/ 



fW2M^fti?tLDf$, Dfffi, DfS, DflbfcU mfffiD f JiCt (ft 
f£D f ffi^^f-t^f^-fi-tfn, |!I^oi#^ii^O|6^ftO/J^v^60^D f 
IStL, mrlfiD f mtBufeD f tW^i^iit'n, 
mm<D'hZ^i><DZD f Mi: U iiuSfiD f jHHtlftiBD f HMbO«S«:^rL 

ft^mz,t U MffiJ^¥^bm^^^*Jci3tt^H(rfS¥MZLO^«¥<7) 

2>M*¥&ft-<0 3 m<Dtp<bX<D f <7)tft£j@;i&^J: ^ CD flEl^'D 

fl£ ¥ <r> itMJj&i <o^k=f- * K is it & ¥ffi 2, (TMi^lS ¥ <£> ^^il^ffl i 19 /h S v ti 
, luf£¥M2 J ^^bm^f6]O^*-ec0'ft^\ »*^Mt^l¥»3 

jjao^^-e^ f <7)m*T®b%^£ 7 ciufeD f mMS.TOiBD f mit* 1 jsit 

[0 0 3 73 

*»WO#=0«HI«l -7°n ^7AT'li, 3 ^^«*±OflUS $ fitz 3 ^7C«it 

1 , m3 ^m&tt<D&m*'bKis^xm s fct -c^a^ ti& 

t coiuisKm f o#iiiiM5Ht<7> 2 m^mz^fi^fiD f Dfs, d f 
Df^b. Df^, Dfjfetu mmv ij&tmmi) fmtmw^m&ii-e 
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ittt, mud f jitmmD f mm^^m^t^v. mw^-^m^^^^m 

fi^Mt^i¥» 7!o^tw f f ^ii^v^J; 9 i:, ItrfSD f jKSu 
buISD f i^b&^fffED f 1 J^T<7)ft^:o^i;J; 0 ^g^t" & £ t tz «t 

Sr*^Tl-&lfi¥W0 7lia^)^t?Of ttfcTEIkfcv*.!:? KU1T8BD f m®> fu 
fr<7>^,£. *m»2f|nJOte^j&, ttM^ft<^T-*u itM^fro^*, 

[0 0 3 83 

[0 0 3 93 

[H l £*fc^JBI] 

[0 0 4 0] 

a f(x, t) / a t =ac(x, t) f(x, t) / a x 
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- (18) 

t -r & mu^-nm^ * m&it l x $m * n -? r°m * # x. %> 0 f u, t ) * stint t £ 

[0 0 4 1 ] 

H[I,J](x) = F(I,J) + G(I.J)x 

- (19) 

&£^£o *»Ww*MW»H[I.J](x)^SS:G(I,J)W^**^MLr*i^"C 
&^& 0 f^-CMJ:^:, G(I, D^^W^HiOi^T^^'J XA<7) 

[0 0 4 23 

s-r & tc, i >xjtmmm±(Dmm $ =m£ i ^7c«ds^ mi xxM&ft* 

f(x,t) t**ffl*$*L/it &m&f(x,t) t&fcf-nmttmm't&ztizxv 
xmm%kF&,T)z, 1 1 ^ii^Hti, j] (x) t Lt^it^fc^^o -t 

[0 0 4 33 
*> L&M&C(I+l/2J)j&*:E^£>*Uf, 

Q(I+1/2,J) := C(I+l/2, J)H[I,J](Ax/2 - C(IJ) At/2) 

- (2 0) 

Q(I+1/2,J) := CU+l/2,J)H[I + lJ](-Ax/2 - C(IJ) At/2) 

- (2 1) 

F(I,J+l) =F(I,J) + At (Q(I+1/2,J) - Q(I-1/2J) )/A x 
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- (2 2) 

[0 0 4 4] 

nfc£, e^hkj i >xmmz}3^xm?.zm-fx u - -fmma, j) ^mro «t ^ k 

^#>&o if, 

Sp( I, J): =F(I+1J) F(I,J) 
Sm( I, J): =F(I,J) F(I-1,J) 

- (2 3) 

AF(I,J) = (sgn(Sp(I,J) + sgn(Sm(I,J)) ) Min(|Sp(I, J) I, |Sm(I, J) I ) 

••• (2 4) 

£SA-f£o {IU sgn () (i^Km-e^^o ^<D? \m&JE<D £ i 

tcti i , ftoBf tcti- i *&-tm%c?:3bz> 0 ^m<Dm[xmom\mmm.^ •? 2 /itt 

^c§v^i:^m^-e^^>o -OAF (I,J)£fiJfflLT 
G(I,J) = AF(IJ)/Ax 

- (2 5) 

[0 0 4 53 

Sp(I,J)/A x(i, iilS^ (2 3) (DTzmfrhMmmfr, Sm(I,J)/Axti±fa 

^ (23) <D7Bmfrb&mmfr$:M.mirz> 0 ^tthtmw^cot^te, (sgn(s P 

(I, J) + sgn(Sm(I,J)) tf^futftlQ, £#7&Mf n £ & £ 0 -e Xte^t £ (i, G 
(IJ)(i. (sgn(Sp(I,J) + sgn(Sm(I,J)) ^2 £ tzit- 2 tli^ 
* * v> (i hu mM53-0^<tO/h £ & 0<z> 2 t —Wet Z> o 
[0 0 4 63 

SHBSfcatr, 4>^T'(7)^^^gi-^^T-o4 I ^-e(7)1itF(X > T)^-^^, l^M 
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JO 0 4 7] 

o 

[0 0 4 83 

L^M^ra^M^F (•, i) zmmit-r&o 

[0 0 4 93 

»f^s2m ±&<dq (-,]) mm^fro (20) rv? 

S (21))). COBf, _hlS^: (2 0) XC^S: (2 1) H^-T J: -9 

fcm<nmmm%tiizmmirz> 0 fit, *7^ys3-cQ (-,j) ^r-*-?^ 
f ( • ,j+n <7)it»^tf^ (±ie^: (2 2) ) o 

[0 0 5 03 

i£±.tm (is) omr B mi^t*i>^ M#^B#p«S7m£ot^?& 

<7)'lfWc|£v^ B#£ijT = A t J fOF^'filF (-,J) F ( • , J -1) j&> <b # 
^7^1^13, B#K*T = A t * (J-l) J&»t>T = A t * JKggHfcff^o 
[0 0 5 13 

fit, '7 7°S 5-Cft^x-r-y ^S^OfW *l*-t&&tnf£\ 
[0 0 5 23 

tcsv^m^t^tii-o -e (D^^m^Ltz <Dtm2 x$>& 0 vt^m 
1 0 0 0 * t- y yi^m-k^m^ uv^ 0 m 2 x^^%wm^wm\n^ 

%zli, m7 -m9l,Z7rsLfzi>CDI<z Jt^ X Wig ftmm ^ftotv^ 
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[0 0 5 33 

«o««B|-c* «9 , 111 4 *±*09S«om*CW«j££^1-:rn 7^it^ 0 
[0 0 5 4] 

i3Mlll4i:^tj:^:, iOfHWSItt, M^CPUl 16^t'Jll 
8f^?>^^ff^i^ A-KfV^^ 1 1 7^7D7k- (S^iSH) r-V 
X^l 1 O^fcfc&fBlt^lSfflk K 1 1 5 ^-*"7X 1 1 4^»bfc&A 

[0 0 5 5] 

i<7)ft^(i> ^t'J 1 1 8_hKH li:^U7n-ft - Mi7&oTl£ftL 

IB*3-£fc>K f^7^^ 1 1 2_tK^$-£& 0 
[0 0 5 63 

[0 0 5 73 
[0 0 5 83 

a F(x,y,t) / 3 t =3 Cx(x,y,t) F(x,y,t) / 8x 
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+ 3 Cy(x,y,t) F(x,y,t) / 9y 

- (2 6) 

t -t h 2 ^tc^ -e oii^^i^-e^> «9 > - n- £ £ mid-it txtmzfto mm 

£#x.£o F(x,y,t) ZmmtL fzi><Dt L-CF(X,Y,T)£#;i&o F^ltX, Y, 

t ti , &TffiWLM & & m&^Jk if mam m * mm-t z>i><?)tirz> 0 

[0 0 5 93 

H[I,J,T](x,y) =F(I,J,T) + Gx(I,J,T)x + Gy(I,J,T)y 

- (2 7) 

&£*E^£o *HM^£^MP^MmH[I, J, T](x,y)(7)#^Gx(I, J, T)&m;y(I,J,T 
)^*a6^^ML-C{im-e^-<^o ^--C&^&idK, Gx(I,J,T)£<7Gy(I,J,T 

) (Dx^-^m^mm^momm r * *v XA^fit^^o g x (i, j, d wgv u 

[0 0 6 03 

2^7n^^t±<7)«^:^^^2^7E:«it1 : »^> ^2^7C«5t^ 

mm-th z t iz x ^xmmma, y,t) i£2 ^x^^±^h[i, j,t] ( x , y ) ^ 

[0 0 6 1 3 

P^tfcCx(I+l/2,J,T)rtMET&tU2\ 
Qx(I+l/2,J,T) := Cx(I+l/2,J,T) 

H[I,J,T](Ax/2 - Cx(I,J,T) At/2, - CyU.J.T) At/2) 

- (2 8) 

Qx(I+l/2,J,T) := Cx(I+l/2,J,T) 

H[I+l,J,T](-Ax/2 + Cx(IJ.T) At/2, - CyU.J.T) At/2) 
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- (2 9) 

Sfc, M^Cy(I,J+l/2,T)^iE-e*ti{i\ 
07(1,1+1/2,1,1) := Cy(I,J+l/2,T) 

H[I,J,T]( - Cx(I,J,T) At/2, Ay/2- Cy(I,J,T) At/2) 

••• (3 0) 

Qy(I,J+l/2,J,T) := Cy(I, J+1/2.T) 

H[I,J+1,T]( - Cx(I,J,T) At/2, 

-Ay/2- Cy(I,J,T) At/2) 

- (3 1) 

F(I,J,T+1) =F(I,J,T) + At (Qx(I+l/2,J,T) - Qx(I-l/2, J,T) )/A x 

+ At (Qy(I,J+l/2,T) - Qy(I,J-l/2,T) )/A y 

- (3 2) 

[0 0 6 2] 
[0 0 6 33 

Sxp( I, J,T): =F(I+1,J,T) — F(I,J,T) 
Sxm( I, J,T): -F(IJ.T) — F(I-IJ.T) 
Syp( I, J,T): =F(I,J+1,T)— F(I,J,T) 
Sym( I, J,T): -F(I,J,T) — F(I.J-l.T) 

- (3 3) 

Sx(I,J) = (sgn(Sxp(I,J,T) + sgn(Sxm(I,J,T)) ) Min( ISxp(I, J, T) I , 
|Sxm(I,J,T)| ) 
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Sy(IJ) = (sgn(Syp(IJ,T) + sgn(Sym(IJ,T)) ) Min(ISyp(I, J,T) I, 
|Sym(I,J,T)| ) 

••• (3 4) 

[0 0 6 4] 

Vppmin = Min(F(IJ+l, T)-F(I,J, T), F(I+1,J+1, T)-F(I,J, T), 
-c F(I+1,J, T)-F(IJ, T) ) 

Vmpmin = Min(F(I-l,J+l, T)-F(I,J, T), F(I,J+1. T)-F(I,J, T), 

F(I-1J, T)-F(I, J, T), -e) 
Vpmmin = Min( -e, F(I+1,J, T)-F(I,J, T) 

F(I,J-1, T)-F(I,J, T), F(I+1,J-1, T)-F(I.J, T)) 
Vmmmin = Min(F(I-l,J, T)-F(I, J, T), -e 

F(I-1,J-1, T)-F(I,J, T), F(I,J-1, T)-F(IJ, T)) 

••• (3 5) 

t, 

Vppmax = Max(F(I,J+l, T)-F(I,J, T), F(I+1,J+1. T)-F(I,J, T), 
-e, F(I+1,J. T)-F(I,J, T) ) 

Vmpmax = Max(F(I-l, J+l, T)-F(IJ, T), F(I,J+1, T)-F(I,J, T), 

F(I-1,J, T)-F(I, J, T), -e) 
Vpmmax = Max( - e , F(I+1,J, T)-F(I,J, T) 

FUJ-I, T)-F(I,J, T), F(I+1J-1 ( T)-F(I.J, T)) 
Vmmmax = Max(F(I-l,J, T)-F(I, J, T), - e 
F(I-1,J-1, T)-F(IJ, T), F(I,J-1, T)-F(I,J, T)) 

- (3 6) 

[0 0 6 51 

hpp(IJ,T) = Sx(I, J, T)/2. + Sx(I, J, T)/2. 
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hmp(IJ,T) = - Sx(I, J, T)/2. + Sx(I, J, T)/2. 
hpm(I,J,T) = Sx(I, J, T)/2. - Sx(I, J, T)/2. 
hmm(I,J,T) = - Sx(I, J, T)/2. - Sx(I, J, T)/2. 

- (3 7) 

[0 0 6 6] 

rphipp=phimp=phipm=phimm=l t LT.fc£^ i> L (hpp > Vppmax ) & 
£>(£\ phipp = VppmaxAipp^tti iti"o i> L^ (hpp < Vppmin ) &£>{i\ ph 
ipp = Vppmin/hppi:S^ l£i"o 
[0 0 6 7] 

4> L (hpm > Vpmmax ) phipm = Vpmmax/hpm <h fi ^ W.~$~o 

3.fzi> (hpm < Vpmmin ) phipm = Vpmmin/hpmh S W.~f~o i> L 

(hmp > Vmpmax ) &£>(£\ phimp = Vmpmax/hmp <bfi 1* it"l"o 
^.fzii L> (hmp < Vmpmin ) &£>ff, phimp = Vmpmin/hmp t S # i£'t"o & L 
(hmm > Vmmmax ) & (£\ phi mm = Vmmmax/hmm <h fi # DHi~ 0 i & L > (hmm 
< Vmmmin ) &£>(i\ phimm = Vmmm i n/hmm t fi § iSi"o J t ^9 ft^Srfr "9 o 
[0 0 6 83 

C<7) «£ -7 LT?#<b tlfcphipp, phimp, phipm, phimmtC^j- LT, 

Gx(I,J, T) = Min (phipp, phipm, phimp, phimm) Sx(I,J,T)/Ax 
Gy(I,J, T) = Min(phipp, phipm, phimp, phimm ) Sy(I,J,T)/Ay 

- (3 8) 

tL/zmz, (i, j) &*\*i<Dmmmm(i,j,T)(x,y)z 

H(I,J,T)(x,y) = F(I,J, T) + Gx(I,J, T) x + Gy(I,J, T)y 

••• (3 9) 

[0 0 6 93 



mtiE#2 003-3067030 



2002-223145 



^-v: 29/ 



[0 0 7 03 

o 

[0 0 7 13 

£Lfc&fcM#ol0»!B*fcF (•, 1) *®mit1-&o 
[0 0 7 23 

^f-yysn-eii, ±jeoq (•, •, j) mmzftd (±is^: (2 

8)^(29)^(30)^(31)), .lOflf, ±fB5£ (28) £<7;£ 

(29) n^ti^:, ^m^mmm^u^mmir^o zlx, xt^si 
3-cq (•, •, j) f (-, • , j + i) mmzfto (_t 

IB^: (22) ) o 
[0 0 7 33 
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m<7)mmz^\ b#£ot = a t j wFoiF (•, j) f (•, 

J - 1 ) ^^T^U, B#F ( g£T= A t * ( J - 1 ) 7^T = A t * J «H 
[0 0 7 4] 

, xr v7°s i 2 ci*), ft^xu/ -7m^^$-ett*^tf do 

[0 0 7 5] 

t^Sv^^fr#^^t^§^o ^IS^£^L£c7)«l6 (a) , (b) 

x^ho 'titmimmm^Jbzmi o (a) , (b) n*tiSL^rai:##-eit 

It^ToT&fK ItfJMi: 1 0 0 0^f7 ylk^m^m^LX^^o El 6 (b 
) ^i-^W^^#rLlt^:^(i, Hi 0 (b) ^LfckOKJfc^-c^lit 
v * - -7 £ ^otw^:^ rt^* & o 

[0 0 7 63 

[0 0 7 73 

m 3 mmrm] 

m 3 Illifli, .LIB* 2 Hit^lg-eKW Lfcr&TCtfM&^fctt-t&WHIT 

8 F(x,y, z, t) / d t =d Cx(x,y, z, t) F(x, y, z, t) / 3x 

+ 8 Cy(x,y, z,t) F(x,y,z, t) / 3y 
+ d Cz(x,y, 2, t) F(x, y, z, t) / 8z 

- (4 0) 
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x.*,F(x,y ( z ( t)!:iMLfc^HtF(X,Y,Z f T)m-6o F<D?mx, Y, 
Z, Tti, rat«^^1MWlitf^t*t-?)^i:t^o 
[0 0 7 8] 

H[I,J,K,T](x,y,z) = F(I,J,K,T) + Gx(I, J.K.T) x 

+ Gy(I,J,K,T)y + Gz(I,J,K,T)z 

••• (4 1) 

mzm^&o ^m^m<ommmm[i,],K,TUx,y,z)<D%mx(i,j,K,T), c y (i 

,J,K,T)SW f Gz(I,J,K,T)0*i6^riC|gLT(±^-eJ6^*o ^-i-CJK^* J: fc, 
Gx(I,J,K,T), Gy(IJ,K,T)^^Gz(IJ,K,T)(D^^^||MSS<7)Mrar^ 
*V XA<n%mx*>Z>o &c(I,J,T)X0 ! Gy(I,J f T)OftseHBRL, »ro«Mir^w 

o 

[0 0 7 93 

o#Jfe^i> is v > -cffi^sE ^ -S^tc iot £ & & & B#£iJ T 3 ^tc« 
ii»^-±OM»F(X, Y, Z, T) 1 sFffiift 3 *tfc £ § , SIW»F(X, Y, Z, T) t Stf»<7>il 

m*mm-r* -^a^ ttw(x, y, z, t) ^k^i • mm txm 3 #sE2iffi«± 

OMtH[IJ > KJ](x ( y,z)^itl»^(7)^^ 0 
[0 0 8 03 
& L&BS IfcCx ( 1+1/2, J , K, T) tfiE X & fUl , 
Qx(I+l/2,J,K,T) := Cx(I+l/2J,K,T) 

H[I,J,K,T](Ax/2 - Cx(I,J,K,T) At/2, 
-Cy(I,J,K,T) At/2,- Cz(I,J,K,T) At/2) 

••• (4 2) 

Qx(I+l/2,J,K,T) := Cx(I+l/2,J,K,T) 

H[I+l,J,K,T](-Ax/2 + Cx(I,J,K,T) At/2, 
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-Cy(I,J,K,T) At/2,- Cz(I,J,K,T) At/2) 

••• (4 3) 

ttz, i> L(>M^cCy(I,J+l/2,K,T)^iE-e*)tt(f, 
Qy(I,J+l/2,K,J,T) := Cy(I.J+l/2,K,T) 

H[I,J,K,T]( - Cx(I,J,'K,T) At/2, 

Ay/2- Cy(I,J,K,T) At/2, - Cz(I, J,K, T) At/2) 

... ( 4 4) 

Qy(I,J+l/2,K,J,T) := Cy(I, J+1/2.K.T) 

H[I,J+1,K,T]( - Cx(I,J,K,T) At/2, 

-Ay/2- Cy(I,J,K,T) A t/2, - Cz(I, J, K, T) At/2) 

— (4 5) 

ifz, i> Li) m$cCz(I, J, K+l/2, T) **jE"C*tL^ 
Qz(I,J,K+l/2,J,T) := Cz(I,J,K+l/2,T) 

H[I,J,K,T]( - Cx(I,J,K,T) At/2, 

- Cy(I,J,K,T) At/2, Az/2- Cz(I,J,K,T) At/2) 

— (4 6) 

Qz(I,J,K+l/2,J,T) := Cz(I,J, K+l/2, T) 

H[I,J,K+1,T]( - Cx(I,J,K,T) At/2, 
-Cy(I,J,K,T) A t/2, -Az/2- Cz(I,J,K,T) At/2) 

••• (4 7) 

F(I,J,K,T+1) =F(I,J,K,T) + At (Qx(I+l/2, J,T) 

-Qx(I-l/2,J,T) )/Ax + At (Qy(I,J+l/2,K,T) 
-Qy(I,J-l/2,K,T) )/Ay+ At (Qz(I, J, K+l/2, T) 
-Qz(I,J,K-l/2,T) )/Az 
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- (4 8) 

[0 0 8 1 3 
[0 0 8 23 

i-t\ 

Sxp( I,J,K,T): =F(I+1,J,K,T) — F(I.J,K,T) 
Sxm( I,J,K,T): =F(I,J,K,T) — F(I-1J,K,T) 
Syp( I,J,K,T): =F(IJ+1,K,T) — F(I,J,K,T) 
Sym( I,J,K,T): =F(IJ,K,T) — F(I,J-1,K,T) 
Szp( IJ.K.T): «F(I.J,K+1,T) — F(I.J,K,T) 
Szm( I.J.K.T): =F(I,J,K,T) — F(I,J,K-1,T) 

— (4 9) 

Sx(I, J,K) = (sgn(Sxp(I, J,K,T)+sgn(Sxm(I, J,K, T)))Min(|Sxp(I, J,K, T) I , 
|Sxm(I,J,K,T)|) 

Sy(I,J,K) = (sgn(Syp(I,J,K,T)+sgn(Sym(I,J,K,T)))Min(|Syp(I,J,K,T)|, 
ISym(I,J,K,T)|) 

Sz(IJ,K) = (sgn(Szp(IJ,K,T)+sgn(Szm(IJ,K,T)))Min(ISzp(I > J f K,T)|, 
ISzm(I,J.K,T)l) 

— (5 0) 

[0 0 8 33 

Vpppmin = Min(F(I,J+l,K, T)-F(I,J,K, T), 
F(I+1,J+1,K, T)-F(I,J,K,T), 

F(I+1,J.K, T)-F(I,J,K, T), 

F(I,J+1,K+1, T)-F(I,J,K, T), 
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F(I+1,J+1,K+1, T)-F(I,J,K, T), 

F(I,J,K+1, T)-F(I,J,K, T),F(I+1,J,K+1, T)-F(I,J,K, T) ) 
Vmppmin = Min(F(I-l, J+l, K, T)-F(I,J,K, T), 
F(I,J+1,K, T)-F(I, J,K, T), 
F(I-1,J,K, T)-F(I,J,K, T), -s, 
F(I-1J+1,K+1, T)-F(I,J,K, T), 
F(I,J+1,K+1. T)-F(I,J,K, T), 
F(I-1,J,K+1, T)-F(I,J,K, T), 
F(I,J, K+l. T)-F(IJ,K, T)) 
Vpmpmin = Min( - e , F(I+1,J,K, T)-F(I,J,K, T) 

F(I,J-1,K, T)-F(I,J,K, T), F(I+1,RKT)-F(I,J,K, T), 
F(I,J. K+l, T)-F(I.J,K, T),F(I+1,J,K+1, T)-F(I,J,K, T) 
F(I,J-1,K+1.T)-F(I,J,K. T), F(I+1,J-1,K+1, T)-F(I, J.K.T)) 
Vmmpmin = Min(F(I-l, J, K, T)-F(IJ,K, T), -e 

F(I-1,J-1,K, T)-F(I,J,K, T), F(I,J-1,K, T)-F(IJ,K, T), 
F(I-1,J,K+1, T)-F(I,J,K, T), F(IJ, K+l, T) -F(I, J, K, T) , 
F(I-1,J-1,K+1, T)-F(I,J,K,T),F(I,J-1,K+1,T)-F(I,J,K,T)) 
Vppmmin = Min(F(I, J+1,K,T)-F(I, J,K, T) , F(I+1, J+l, K, T)-F(I, J,K,T) , 

F(I+1,J,K, T)-F(I,J,K, T), 
F(I,J+1,K-1,T)-F(I,J,K, T),F(I+1,J+1,K-1, T)-F(I, J,K,T), 
F(I,J,K-1, T)-F(I,J,K, T), F(I+1,J,K-1, T) -F(I, J, K, T) ) 
Vmpmmin = Min(F(I-l, J+l, K, T)-F(I, J, K, T) , F(IJ+1, K, T)-F(I, J,K,T) , 
F(I-1,J,K, T)-F(I,J,K, T), -e, 

F(I-1,J+1,K-1,T)-F(IJ,K,T), F(I,J+1,K-1, T)-F(I, J.K.T), 
F(I-1J,K-1, T)-F(I,J,K, T), F(I,J, K-l, T)-F(I, J,K,T)) 
Vpmmmin = Min( -e, F(I+1,J,K, T)-F(I,J,K, T) 

F(I,J-1,K, T)-F(I,J,K, T), F(I+1,J-1,K T)-F(I,J,K, T), 
F(IJ,K-1, T)-F(I,J,K, T), F(I+1,J,K-1, T)-F(IJ f K, T) 
F(IJ-1,K-1,T)-F(I,J,K,T), F(I+1.J-1,K-1, T)-F(I, J.K.T)) 
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Vmmmmin = Min(F(I-l,J,K, T)-F(I,J,K, T), -e 

F(I-1,J-1,K. T)-F(I,J,K, T), F(I,J-1,K, T) -F(I, J,K, T), 
F(I-1,J,K-1, T)-F(I,J,K, T), F(I,J, K-l, T)-F(I,J,K, T), 
F(I-1.J-1,K-1.T)-F(I,J,K,T), F(I,J-1,K-1. T)-F(I, J,K,T)) 

- (5 1) 

Vpppmax = Max(F(IJ+l,K,T)-F(I,J,K,T), F(I+1,J+1,K, T)-F(I, J.K.T), 

F(I+1,J,K, T)-F(IJ,K, T). 
F(I,J+1,K+1,T)-F(I,J,K,T),F(I+1,J+1.K+1,T)-F(I.J,K,T), 
F(I,J,K+1, T)-F(I,J,K,T),F(I+1.J,K+1.T)-F(I,J,K,T)) 
Vmppmax = Max(F(I-l, J+1,K,T)-F(I, J,K,T),F(I, J+l.K, T)-F(I, J,K,T), 
F(I-1,J,K, T)-F(I,J,K, T), - e , 

F(I-1,J+1,K+1,T)-F(I,J,K.T),F(I,J+1,K+1, T)-F(I, J.K.T) , 
F(I-1J,K+1,T)-F(I,J,K,T),F(I,J,K+1, T)-F(IJ,K, T)) 
Vpmpmax = Max( - e , F(I+1,J,K, T)-F(I,J,K, T) 

F(I,J-1,K, T)-F(I,J,K, T), F(I+1,J-1.K T)-F(I,J,K, T), 
FQJ.K+1, T)-F(I,J,K, T),F(I+1.J,K+1, T)-F(I,J,K, T) 
F(I,J-1.K+1,T)-F(I,J,K, T),F(I+l,J-l f K+l,T)-F(I f J.K.T)) 
Vmmpmax = Max(F(I-l, J, K, T)-F(I,J,K, T), -e 

F(I-1,J-1,K, T)-F(I,J,K, T). F(I,J-1,K, T)-F(I,J,K, T), 
F(I-1,J,K+1, T)-F(I,J,K, T), F(I,J, K+l. T)-F(I, J.K.T), 
F(I-1.J-1.K+1,T)-F(I,J,K.T).F(I,J-1,K+1, T)-F(I, J.K.T)) 
Vppmmax = Max(F(I, J+l.K, T)-F(I, J,K, T), F(I+1, J+l, K, T) -F(I, J, K, T), 

F(I+1J,K, T)-F(I,J,K, T), 
F(I,J+1,K-1, T)-F(I,J,K,T),F(I+1,J+1,K-1, T)-F(I, J.K.T), 
F(I,J,K-1, T)-F(I,J.K, T), F(I+1,J,K-1, T)-F(I, J,K,T)) 
Vmpmmax = Max(F(I-l, J+1,K,T)-F(I, J,K,T),F(I, J+l.K, T)-F(IJ,K, T), 
F(I-1,J,K, T)-F(I,J,K, T), - e , 

F(I-l f J+l,K-l,T)-F(I,J,K,T),F(IJ+l,K-l > T)-F(I, J.K.T), 
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F(I-1,J,K-1,T)-F(IJ,K,T),F(I,J,K-1,T)-F(I,J,K, T)) 
Vpmmmax = Max ( - e , F(I+1,J,K, T)-F(I,J,K, T) 

F(I,J-l f K, T)-F(I,J,K, T), F(I+1,J-1.K T)-F(IJ,K, T), 
F(I,J,K-1, T)-F(I.J,K, T), F(I+1, J.K-l, T)-F(I, J,K, T) 
F(I,J-1,K-1,T)-F(IJ,K,T),F(I+1,J-1,K-1, T)-F(I. J.K.T)) 
Vmmmmax = Max(F(I-l, J, K, T)-F(I,J,K, T), -e 

Fd-l.J-l.K, T)-Fd,J,K. T), Fd.J-l.K, T)-F(I,J,K, T), 
Fd-1,J,K-1, T)-Fd,J,K, T), Fd,J, K-l, T)-F(I,J,K, T), 
Fd-l.J-l.K-l,T)-Fd,J.K,T),F(I,J-l,K-l, T)-Fd,J,K, T)) 

••• (5 2) 

[0 0 8 4J 

















hppp(I,J,K,T) = 


Sxd,J,K, 


T)/2. + 


Sxd.J.K, 


T)/2 h 


- Sz(I,J,K,T)/2. 


hmpp(I,J,K,T) = 


- Sxd.J.K. 


T)/2. + 


Sxd.J.K, 


T)/2 


+ Sz(I,J,K, 


T)/2. 


hpmpd.JXT) = 


Sxd.J.K, 


T)/2. - 


Sxd.J.K. 


T)/2 


+ Sz(I,J,K, 


T)/2. 


hmmp(I,J,K,T) = 


- Sxd.J.K, 


T)/2. - 


Sxd.J.K, 


T)/2 


+ Sz(I,J,K, 


T)/2. 


hppmd,J,K,T) = 


Sxd.J.K, 


T)/2. + 


Sxd.J.K, 


T)/2 


- Szd.J.K, 


T)/2. 


hmpm(I, J,K,T) = 


- Sxd.J.K, 


T)/2. + 


Sxd.J.K, 


T)/2 


- Szd.J.K, 


T)/2. 


hpmm(IJ,K,T) = 


Sxd.J.K, T)/2. - 


Sxd.J.K, 


T)/2 


- Sz(I,J,K, 


T)/2. 


hmmmd, J,K,T) = 


- Sxd.J.K, 


T)/2. - 


Sxd.J.K, 


T)/2 


- SzdJ.K, 


T)/2. 



- (5 3) 



[0 0 8 5] 

H K > fph i ppp = ph i mpp=ph i pmp=ph i mmp= ph i ppm = ph i mpm=ph i pmm=ph i mmm= 1 b 
LTj3^ > & L (hppp > Vpppmax ) &£>t£\ phi ppp = Vpppmax/hppp b S iff 
"t"o i> L (hppp < Vpppmin ) &£>{£\ phippp = Vpppmin/hppptS^ i£i"o 

& L (hpmp > Vpmpmax ) &£>(3\ phipmp = Vpmpmax/hpmp b S § fii"o 
£ (hpm < Vpmmin ) &£>(£\ phipm = Vpmpmin/hpm<b fit # iSi"o & L 
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(hmpp > Vmppmax ) &<b(i\ phimpp = Vmppmax/hmpp<hg^ilLi~o i tz i> L> 
(hmpp < Vmppmin ) &£>(£\ phimpp = Vmppmin/lmpp <h fi ^ SE^o & L (hmmp 
> Vmmpmax ) phimmp = Vmmpmax/hmmp<hS§Ifi^o 

[0 0 8 6] 

i & L> (hmmp < Vmmpmin ) phimmp = Vmmpmin/hmmp<hflr: 

o & L (hppm > Vppmmax ) phippm = Vppmmax/hppm <h il § ifii^o £ £ 

i> L (hppm < Vppmmin ) &£>(i\ phippm = Vppmmin/hppm<h fli § iSi~o & L ( 
hpmm > Vpmmmax ) phipmm = Vpmmmax/hpmm hl^Sto ^ fzi> (h 

pmm < Vpmmmin ) phipmm = Vpmmmin/hpmmt fi ^ iSi^o & L (hmpm > 

Vmpmmax ) phimpp = Vmppmax/hmppchg^iiE^o 

[0 0 8 73 

^ 7t 4> L> (hmpm < Vmpmmin ) phimpm = Vmpmmin/hmpm<t WL ^ Hi" 

o & L (hmmm > Vmmmmax ) & ff\ phimmm = Vmmmmax/hmmm <b g ^ ifii~ 0 itz 
i> L> (hmmm < Vmmmmin ) phimmm = Vmmmmin/hmmm<h fi ^ itti"o J t 

[0 0 8 83 

j: ^ I: LT|#^> tifcphippp, phimpp, phipmp, phimmp, phippm, phimpm, phi pm 
m,phimmmtc£t LT, 

Gx(I, J, K, T) =Min(phippp, phimpp, phipmp, phimmp, phippm, phimpm, phipmm, phimm 
m)Sx(I,J,K,T)/Ax 

Gy (I, J, K, T) =Min(phippp, phimpp, phipmp, phimmp, phippm, phimpm, phipmm, phimm 
m)Sy(I,J,K,T)/Ay 

Gz (I, J, K, T) =Min(phippp, phimpp, phipmp, phimmp, phippm, phimpm, phipmm, phimm 
m)Sz(I,J,K, T)/Az 

t Ltz'ikK, (I,J,K) ^0«BH»H(I,LKJ)(x,y,z)S: 
H[I,J,K,T] (x,y,z) 

= F(I,J,K,T) + Gx(IJ,K,T)x 
+ Gy(I,J,K,T)y + Gz(I,J,K,T)z - (5 4) 
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[0 0 8 9] 

W2!5Mt0 2f&0|t£**t^:ft Dfl, Dfl, Dfi, Dft, Df^ D 
-tfn, Plff Ot^lif ^i<7)/h$i40^D f ^it U DffltDf 

*^ri"*fl&¥^^7fl&ott»^-cw f flt^^x.^v^9 k> d f$HMD f m 

[0 0 9 0] 

nLxfs^tzmt*. imttm<D&j-&, mm^mco^^ mMttfov 

^bS%^[^o^T-*> mM^ft<7>t&T*U 2fctf 

btitz¥m*mm¥ftx<Dm$mmmmt Lr^^m^x t)mmm%t*&z> 

o ^ co r ;v =r u x a t Jtie^ t £ r ;u n'>; x^t mffi-c& z> 0 

[0 0 9 1 ] 

i>(Dtmmx&z> 0 ttz, *n3mm%m<vm$iim*m3Jkvm4xmmLtz 
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[0 0 9 21 

^ f> & ^ * ± # x. e> ft <b 2 ^tc h h v > l± 3 %Jk& m ± O M & £ 

^!lx.fflftT-*-^<7)S^, 2 ^jc*4 v>li 3 ^x<0f v ^itt©f - ^ <0 
[0 0 9 33 

^tiCPU^MPU) ^f2«^tcf^$^fc-7°n^9A3- K*tt*WL2Sfr 
[0 0 9 4] 

m^mm z mm.1- h £ t k * 0 , *o ■/ n ^ a 3 - k * i sn t £ 

(i, #!jxJ£\ 7n 7 e- (g&iSfcQ T4*9^ ^-FfU^ $17*4*?, 
jm.%7*4*C CD-ROM, CD-R, i&$C^-7\ ^JfHteO* * V # - 

^mt^^M £ ft 2> i> & * ft & - t iim -9 £ -C & v > c 
[0 0 9 53 
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*v\zm&&z*itz'^ wd^7A3- Kojg^t^o^, zvntMmt&t: 

[0 0 9 63 

(V0F&) Kfcv*T, mR*i'*-'7l l zfkW2<£t:tZ1$m&m2-£&Z 

[mi 3 

[0 2 3 
[H3 3 

[H43 

[H53 
[H63 

^2^TO3^Wm^M^^^fUfflL7t:^rLft#0|g*^^llT^^o 
[H7j 

& * <r> m mil & * * "j ffl l fc ^ % it # <r> m * * ^ •* m -c * * „ 

[H8 3 
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